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Abstract. Entanglement in the ground state of the XY model on the infinite chain 
can be measured by the von Neumann entropy of a block of neighboring spins. 
We study a double scaling limit: the size of the block is much larger than 1 but 
much smaller than the length of the whole chain. The entropy of the block has an 
asymptotic limit. We study this limiting entropy as a function of the anisotropy and 
of the magnetic field. We identify its minima at product states and its divergencies 
at the quantum phase transitions. We find that the curves of constant entropy are 
ellipses and hyperbolas and that they all meet at one point (essential critical point). 
Depending on the approach to the essential critical point the entropy can take any 
value between and oo. In the vicinity of this point small changes in the parameters 
cause large change of the entropy. 
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1. Introduction 

Entanglement is a primary resource for quantum information processing [TJ El El Sl- 
it shows how much quantum effects we can use to control one system by another. 
Stable and large scale entanglement is necessary for scalability of quantum computation 
[5J, [6] . The entropy of a subsystem as a measure of entanglement was discovered in [3]. 
Essential progress has been achieved in the understanding of entanglement in various 
quantum systems El El US [III CGI [131 El ESI EE [13 [IE1 [IS EQl Ell E21 E31 ES ESI 
QUEUES. 



The importance of the XY model for quantum information was emphasized in 



[IH EHl EH E2jil- In this paper we consider the entropy of a block of L neighboring 
spins in the ground state of the XY model [on the infinite chain] in the limit L — ► oo. 
We use the results of [331 EU ES] and extend them to the whole phase diagram of the 
modeljj The Hamiltonian of the XY model is 

oo 

n = - J2 0- + tVX+i + (i - + h < (i) 

n=— oo 

Here 7 > is the anisotropy parameter; a^, o\ and a z n are the Pauli matrices and 
h > is the magnetic field. The model is clearly symmetric under the transformation 
7 — > —7 or h — > — h. In [33l [35] only the case < 7 < 1 was discussed, here we can 
confirm that those results can be directly extended for 7 > 1 by analytical continuation. 

The XY model was solved in [351 EZ1 EH1 ES]- The methods of Toeplitz 
determinants and integrable Fredholm operators were used for the evaluation of 
correlation functions, see [3H HQl HH H21 H3J 03]. The idea to use the determinants 
for the calculation of the entropy was put forward in [To] . 

I It is interesting to note that the critical behavior of the XY model is similar to the Lipkin-Meshkov- 
Glick model [31]. 

§ Moreover, in the Appendix we explain how these results can be used to calculate the entanglement 
of the XY model in a staggered magnetic field. 
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The solution of the XY model looks differently in three cases: 
Case 2 is defined by h > 2: This is strong magnetic field. 

Case 1a is defined by h < 2 and 7 > a/1 — (h/2) 2 : Moderate magnetic field for 
small anisotropy, and includes zero magnetic field for large anisotropy. 

Case 1b is defined by h < 2 and 7 < a/1 — (h/2) 2 : It describes weak magnetic 
field, including zero magnetic field in the small anisotropy regime. 

At 7 = and for h < 2 the model is known as the isotropic XY model (or XX 
model) and its spectrum is gapless. The entropy for this critical phase was calculated 
in [16]. The boundary between cases 2 and 1a (h = 2) is also critical. In the rest phase 
diagram, the spectrum of the XY model is given by: 

e k = \J (cos k - h/2) 2 + 7 2 sin 2 k. (2) 

We draw the phase diagram and the three cases we are considering in Fig. [IJ 

At the boundary between cases 1a and 1b (h = 2a/1 — 7 2 ) the ground state can 
be expressed as a product state as it was discovered in [45J . The ground state is in fact 
doubly degenerated: 

lOSiH J] [cos(0)|U>+sin(0)U n >], 
nelattice 

\GS 2 ) = H [cos(0)| Tn) - sin(0)| U] (3) 
nelattice 

Here cos 2 (2#) = (1 — 7)/(l + 7). The role of factorized states was emphasized in 
[40j HU HT] . Let us mention that the rest of the energy levels are separated by a gap 
and correlation functions decay exponentially. The boundary between cases 1a and 1b 
is not a phase transition. 
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Figure 1. Phase diagram of the anisotropic XY model in a constant magnetic field 
(only 7 > and h > shown). The three cases 2, lA, 1b, considered in this paper, 
are clearly marked. The critical phases (7 = 0, h < 2 and h = 2) are drawn in bold 
lines (red, online). The boundary between cases lA and 1b, where the ground state 
is given by two degenerate product states, is shown as a dotted line (blue, online). 
The Ising case (7 = 1) is also indicated, as a dashed line. 
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2. Block Entropy 

In general, we denote the ground state of the model by \GS). We consider the entropy 
of a block of L neighboring spins: it measures the entanglement between the block 
and the rest of the chain [51 [3U]. We treat the whole ground state as a binary 
system \GS) = \A&cB). The block of L neighboring spins is subsystem A and the 
rest of the ground state is subsystem B. The density matrix of the ground state is 
Pab = \GS)(GS\. The density matrix of the block is pa = Ttb(pab)- The entropy 
S(pa) of the block is: 

S(p A ) = -Tr A (p A \np A ) (j (4) 

Note that each of the ground states ((3]) is factorized and has no entropy. 

To express the entropy we need the complete elliptic integral of the first kind, 

I(Jfe) = t d% = * (5) 

Jo v/(l-o; 2 )(l-ft 2 ) 

and the modulus 

T = I(k')/I(k), k' = Vl^W (6) 

The magnetic field and anisotropy define the elliptic parameter k: 
7 / vW 2 ) 2 + 7 2 - 1, Case 2 



k = < v/(V2) 2 + 7 2 -l / 7, Case 1a (7) 

^ y/l -7 2 - (h/2) 2 j y/l - (h/2) 2 , Case 1b 
k vanishes at large magnetic fields (h — > oo), at 7 = for h > 2 and at the boundary 
between case 1a and 1b (h = 2yl — 7 2 ). In all these regions of the phase diagram, 
the ground state of the system is given by product states (a ferromagnetic state in the 
first two cases, the doubly degenerate state for the latter). At the phase transitions 
(h = 2 and 7 = 0, h < 2) the elliptic parameter k = 1. 
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Figure 2. The limiting entropy as a function of the magnetic field at constant 
anisotropy 7 = 1/2. The entropy has a local minimum S = In 2 at h = 2yl — 7 2 
and the absolute minimum for h — ► 00 where it vanishes. S* is singular at the phase 
transition h = 2 where it diverges to +00. The three cases are marked. 



In the paper [53] we used determinant representation for the evaluation of the 
entropy. The zeros of the determinant form an infinite sequence of numbers: 



, , 1-0" 
tanh I m H — | 7tt 



A, 

here a = 1 in Case 1 and a = in Case 2. Note < X m < 1 and A m — >• 1 as m — > 00 

These zeros allowed us to represent the entropy as a convergent series in |33j : 

2 



9 



S(p x )= E (! + A -) ln Y 



m=— 00 



+ A r , 



(9) 
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I. Peshel also obtained the series in cases of non-zero magnetic field, see [34J[|J] He 
summed it up into: 
1 



S(pA) = ~ 

S(pa) = I 



In 



Case 2 (10) 

k> \ (n 2I{k)I{k') 



7T 



+ ln 2, Case 1a (li; 



16k' 

In our paper [33], we have shown that equation (Q is valid for all three cases, which 
allowed us to sum up the series © in case of weak magnetic field (including zero 
magnetic field) as well: 



* S'(pa) = ~ 



In ( + ( 2 - k 2 ) 

\16k' J V y 



+ ln 2, Case 1b, (12) 



7T 

The rigorous proof and the precise history is given in the paper |35J. Here we report 
that these results are valid also for 7 > 1. 

3. Entropy's analysis 

Now we can study the range of variation of the limiting entropy. We find a local 
minimum S(pa) = In 2 at the boundary between cases 1a and IB (h = 2^1 - 7 2 ). 
This is the case of doubly degenerated ground state (j3J) and it is consistent with [35] , 
where it was shown that when the ground state becomes a superposition of two product 
states with different quantum numbers, then the entropy of a subsystem turns into In 2. 

The absolute minimum (S = 0) is achieved at infinite magnetic field or at 7 = 
for h > 2, where the ground state becomes ferromagnetic (i.e. all spins are parallel). 
The entropy diverges to +00, i.e has singularities, at the phase transitions: h = 2 or 
7 = 0. To show this behavior of the limiting entropy, we plot it as a function of the 
magnetic field h at constant anisotropy 7 = 1/2 in Fig. [2J 

|| In comparing with the results of [34] . the reader should keep in mind that Peshel calculated the 
entropy per boundary, therefore his results differ by a factor of 2 compared to the ones in this paper. 
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Figure 3. The limiting entropy as a function of the anisotropy parameter at constant 
vanishing magnetic field h = 0. The entropy has a minimum 5* = In 2 at 7 = 1 
corresponding to the boundary between cases 1a and 1b. S diverges to +00 at the 
phase transition 7 = 0. 



We provide a plot of the entropy as a function of 7 at constant (vanishing) magnetic 
field (h = 0) in Fig. [31 where the local minimum S = In 2 is again visible at 7 = 1. 
Figure H] is a three-dimensional plot of the entropy as a function of both the anisotropy 
parameter 7 and of the external magnetic field h and all the feature discussed so far 
are visible. 

The Ising point. The degenerate product states case ([3]) is particularly 
interesting. For h — 0, 7 = 1, the XY model reduces to the Ising model and the 
ground state is given by the Bell-Pair state: 



ICS,) 



n ^(iu+(-inu) 



nglattice 
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Figure 4. Three-dimensional plot of the limiting entropy as a function of the 
anisotropy parameter 7 and of the external magnetic field h. The local minimum 
S = In 2 at the boundary between cases 1a and 1b is visible and marked by a 
continuum line. S diverges to +00 at the phase transitions h = 2 and 7 = 0, h < 2. 
The entropy takes every positive value in the vicinity of the essential critical point 
(^,7) = (2,0). 




We plot the entropy as a function of the magnetic field at the Ising point 7 = 1 in 
Fig. One can notice that the local minimum S = In 2 is achieved at h = 0. 

The essential critical point. Another interesting limit is reached around the 
point 7 = 0, h = 2. This point belongs to both of the critical phases of the XY model, 
so the entropy does not have an analytical expression [fixed value] on this point, but 
we can study the behavior of S(pa) in vicinity of this point. We already studied a 




k=0 

' 1 ' - h 

12 3 4 

Figure 5. The limiting entropy as a function of the magnetic field at the Ising point 
7=1. The entropy has a local minimum S = In 2 at h — and the absolute minimum 
for h — > oo where it vanishes. S is singular at the phase transition h = 2 where it 
diverges to +oo. 

couple of trajectories reaching this critical point: along the boundary between cases 
1a and 1b (h = 2yl — 7 2 ) the entropy is on its local minimum S = In 2. Along the 
critical lines (h = 2 and 7 = for h < 2) the entropy is divergent, while for 7 = and 
h > 2 the ground state is ferromagnetic and the entropy is 0. Since the limit of the 
entropy reaching the point (h,j) = (2,0) does not exist (it is direction-dependent), we 
call this point essential critical point. In the next section we study the vicinity of this 
point and show that depending on the direction of approach the entropy can take any 
positive value. In Fig. [6] we present a contour plot of the entropy around the essential 
critical point. From this plot one can see that the entropy can assume a wide range of 
values near the point. 
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Figure 6. Contour plot of the limiting entropy near the essential critical point h = 2, 
7 = 0. Regions of similar colors have similar entropy values and the lines where colors 
change are lines of constant entropy. S(pa) is diverging to +oo on the critical lines 
h = 2 and h < 2, 7 = 0. One can see that near the essential critical point the lines of 
constant entropy grow denser. 



4. Ellipses and hyperbolas of constant entropy 



We now look for curves of constant entropy. Since the entropy depends only on the 
elliptic parameter (171) , curves of constant entropy are curves of constant k. Such 
trajectories are easily found and the family of curves of constant entropy can be written 
in terms of a single parameter k: 



Case 2 



Case 1a 



Case 1b 



h > 2 : 
h<2, 

7 > VI " (V2) 2 
h < 2, 



7 < V 1 ~ (h/2f 



DAI 



1, < k < oo (14) 
1, K>1 (15) 

1, K<1. (16) 
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For h > 2 (Case 2), the curves of constant entropy are hyperbolas, while for h < 2 
(Case 1a and 1b) they are ellipses. 

Each point in the phase diagram of the XY model belongs to one of these curves. 
By selecting a value of the parameter k, we select a family of point with the same 
elliptic parameter k in ([7]). There is a one-to-one correspondence between k and k: 



K 2 



k = \lT^ *' = VlT7? Case2 (17) 



1 K 2 - 1 I 

k = \l — — k' = - Case 1a (1* 



K 2 K 



k = Vl-n 2 k' = K Case 1b. (19) 

We recognize that k = 1, as the boundary between Case 1a and 1b, is the curve where 
the ground state can be expressed as a doubly degenerate product state (J3j). 

It is important to notice that this curves have all the essential critical point 
(h, 7) = (2, 0) in common. This means that starting from any point in the phase 
diagram of the XY model, one always reaches the essential critical point by following 
a curve of constant entropy. 

For h < 2, the entropy has a minimum at In 2 and diverges to +00 approaching 
the critical line h = 2. For h > 2, S(pa) decreases monotonically from +00 near the 
critical line to at infinite magnetic field. Beside the critical lines, the entropy is a 
continuous function, so its range is the positive real axis. 

This means that, depending on the direction of approach, the entropy assumes 
every positive number near the essential critical point, since every ellipse or hyperbola 
of constant entropy passes through that point. In other words, a small deviation from 
the essential critical point can bring a big change in the value of the entropy. This 
is very important from the point of view of quantum control, because it allows to 
change dramatically the entanglement (and hence the quantum computing capabilities) 
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with small changes in the parameters of the system. 

It is easy to express the entropy in terms of the parameter k defining the ellipses 
and hyperbolas of constant entropy: 



S(pa) 

S(pa) 
S(pa) 



^4(^+1) 2^ 



In 



In 



K 7T K 2 + 1 

- 2 - 1 2 U 2 + 1) r 



16k 



16k 



+ -- 



TT K 



-(k 2 + 1)I(VT^ )I(k) 

TT 



5. Entropy approaching the critical lines 



K 2 + 1 



II- 

K 



+ ln 2, 



Case 2 



+ ln 2, Case 1a 



Case 1b. 



(20) 



Using the formulae of the previous section, we are now in a position to discuss the 
divergences of the entropy near the critical phases. General results exist in these cases, 
based on a conformal field theory approach [20] and specific results were derived for 
the isotropic case (7 = 0) in [16j. We know that, in the double-scaling limit we 
are considering, the entropy diverges logarithmically with the size of the block. The 
coefficient of this logarithmical divergence can be calculated by knowing the central 
charge of the corresponding conformal field theory at the critical point [20J. 

Setting k = or k = 00, the ellipses and hyperbolas of constant entropy f|T6|) 
collapse into the critical lines, i.e. a vertical line at 7 = or a horizontal line at h = 2, 
respectively. Using (|20|) . we can study how the entropy diverges approaching these 
lines. Using Case 1b in (|20|) . we can take k — > to find 

1 



S(K^0,h< 2) 



i,„(2) + i>„[i-( ft /2)T + ... 



(21) 
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which is consistent with the results obtained in [TB] for the isotropic case (7 = 0). 

We can investigate how the entropy approaches the critical line h = 2 from below 
and from above. In the former case, we shall set k — > 00 in Case 1a of ( |2"U1) : 



In the latter case, for a direction almost parallel to the critical line h = 2, but slightly 
above it, we take k — > 00 in Case 2 of (120]) : 



These results are in agreement with the conclusion of [2D] . 
6. Conclusions 

We analyzed the entanglement in the ground state of the XY model on the infinite 
chain by studying the von Neumann entropy S(pa) of a block A of neighboring spins. 
This entropy is an effective measure of the quantum computing capabilities of a system 
and plays a fundamental role in the field of quantum information. 

Using previously known results for the entropy in the limit of a large block of spins 
[35]-[33], we studied the behavior of S(pa) in the phase diagram of the XY model (see 
Table 1). We found that for h < 2, the entropy has a local minimum S = In 2 on the 
curve (h/2) 2 + 7 2 = 1. On this line the ground state is a doubly degenerate linear 
combination of product states. The entropy diverges to +00 at the phase transitions 
h = 2 and 7 = 0, h < 2. For h > 2, the entropy reaches the absolute minimum 
at infinite magnetic field h — > +00 and for 7 = 0, i.e. when the ground state is 




(22) 



S(k -> 00, h > 2) ~ - In (4 k) + . . . 
3 



[(V2) 2 -l]+- ln(4 7 ) + ... 



(23) 
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Region 




Curves of 
Constant S 


Range of 
Parameters 


2: 


h>2 


1 

6 


[ 1 4 , 2(fc 2 -fc' 2 )/(fc)J(fc')l 

111 k k' "'" 7T 


/ i* \ 2 /-v\2 

(1) -O =1 


< fc < 1 
U ^ K < oo 


H 


f h<2, 

[ 7 > \A " (V2) 2 


l 
6 


fca 2(2-fc 2 )/(fc)/(fc')" 

111 16Jfe' ^ 7T 


+ In 2 


(l) a + G) 8 = i 


< fc < 1 

K > 1 


H 


' /i< 2, 

w 7 < Vl - (V2) 2 


1 
6 


"i fc2 2(2-fc 2 )/(fc)/(fc'f 

1U 16Jfc' " 1 ~ 7T 


+ In 2 


(j) a +e) 2 =i 


< fc < 1 
« < 1 

fc = Vl - 


7 = y/1 ~ {h/2f 


In 2 


(|) 2 + 7 2 = 1 


fc = 
«; = 1 



Table 1. Recap of the results with the entropy in the different regions of the phase 
diagram, the curves (ellipses and hyperbolas) of constant entropy and the relationship 
between the elliptic parameter k and the parameter k defining the family of curves. 



ferromagnetic. S(j>a) diverges to +oo on the critical line h = 2 and it is continuous 
otherwise. 

We identified a set of curves (ellipses and hyperbolas) of constant entropy. They 
are given in (11411161) . All these curves have one point in common, that we decided to 
call essential critical point: (h,j) = (2,0). The fact that all the curves of constant 
entropy pass through one point, together with the fact that the range of the entropy 
as a function of 7 and h is the positive real axis, means that the entropy can assume 
any real positive value near the essential critical point, depending on the direction of 
approach. In turn, this means that the essential critical point is very important for 
quantum control, in that small changes in the parameters can change the entanglement 
dramatically. 
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With this work, we conclude the analysis of the asymptotic Von Neumann entropy 
for the bi-partite one-dimensional XY model. We covered the whole phase diagram 
(including 7 > 1), focusing on the sector h > and 7 > 0: since the model is invariant 
under the substitution 7 — > —7 or h — > —h the results for the entanglement can 
be extended immediately to negative values of the anisotropy parameter 7 or of the 
magnetic field h. 

Finally we note that the work done so far on the XY model in a constant magnetic 
field allows us to calculate the bipartite entropy of the XY model in a staggered 
magnetic field as well. As we discuss in the appendix, there is exact mapping between 
these two models. Therefore, the knowledge of the entanglement for one of the models 
automatically gives the entanglement for the other. We give some details in the 
appendix. 
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Appendix A. XY model in a staggered magnetic field 

It is a well known fact in the theory of integrable models, that there is an exact 
mapping between the traditional XY model in constant magnetic field described by 

the Hamiltonian ([1]) and the XY model in staggered magnetic field: 

00 

W = - J C 1 + 7>X+i + (1 - + (-1)" h' <• (A.l) 
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This mapping is achieved by performing a rotation of every other spin alongthe x 
direction. To identify the two Hamiltonians, one also needs to substitute 7 — ► 1/ 7j]j and 
to rescale the magnetic field and the Hamiltonian by a factor of I/7 and 7 respectively: 

7=1/7 ti = h/~f J = 7. (A.2) 

In the main body of this article we analyzed the entanglement of the XY model 
in a constant magnetic field ([T]). The results we derived can be applied directly to 
calculate the bi-partite entanglement of the XY model in a staggered field (IA.1I) . All 
formulae are valid and to calculate the entropy for a staggered field one only needs to 
take the appropriate result and perform the substitutions (1A.2jl . 

Using (1A.2I) . the spectrum of the XY model in a staggered magnetic field is: 



e k = ^(ycosk-h'^y + sin 2 k. (A.3) 

From this, we see that the critical phase h = 2 is mapped to the line h' = 27'. 

Therefore, the mapping of the different cases and the definitions of the elliptic 
parameter for this model are: 



Case 2 



I ti > 2 7 ' : k = l/y/(h'/2) - 7 ' 2 + 1 (A.4) 



ti < 2 7 ', 



Case 1a <^ ' : k = y/(h'/2) - i 2 + 1 (A.5) 

V< y/l+{h'/2f 

, ti < 2 7 ', /y2 _ (h'/2) 2 - 1 

Case 1b ' J : k = v \ y 1 } =- . (A.6) 



7' > VTTWJW ' V"?' 2 - (h'/2y 

We draw the phase diagram of the XY model in a staggered field and indicate 
the three cases in Fig. IA11 With the definitions (IA.4HA.6l) . one can plug the elliptic 



parameter into (11 OH 121) and use the other results of this paper to calculate the entropy 
of the XY model in a staggered magnetic field. 

Note that the large anisotropy regime is so mapped into the small anisotropy regime and vice- versa. 
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Figure Al. Phase diagram of the anisotropic XY model in a staggered magnetic 
field (only 7' > and h' > shown). The three cases 2, lA, 1b, considered in this 
paper, are clearly marked. The critical phase h! = 27') is drawn as bold line (red, 
online). The boundary between cases lA and 1b, where the ground state is given by 
two degenerate product states, is shown as a dotted line (blue, online). 
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